The origin of the E + transition in GaAsN alloys 
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Optical properties of GaAsN system with nitrogen concentrations in the range of 0.9-3.7% are 
studied by full-potential LAPW method in a supercell approach. The E+ transition is identified 
by calculating the imaginary part of the dielectric function. The evolution of the energy of this 
transition with nitrogen concentration is studied and the origin of this transition is identified by 
analyzing the contributions to the dielectric function from different band combinations. The Lie- 
derived states are shown to play an important role in the formation of the E+ transition, which 
was also suggested by recent experiments. At the same time the nitrogen-induced modification of 
the first conduction band of the host compound are also found to contribute significantly to the E+ 
transition. Further, the study of several model supercells demonstrated the significant influence of 
the nitrogen potential on the optical properties of the GaAsN system. 

PACS numbers: 71.15-m, 71.45.Gm, 71.55.Eq 
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GaAsi_ 2 ,N 2; alloys have attracted much attention of 
both experimentalists and theoreticians due to their un- 
usual physical properties. Among those features, ob- 
served in the typical nitrogen concentration regime of 
0.5%-3.0%, one should mention the anomalous band gap 
reduction, which makes this system technologically im- 
portant for such electronic devices as infrared diode 
lasers 1] and multijunction solar cells Another strik- 
ing feature of this system is the appearance of a new 
composition-dependent optical transition, called E + , de- 
tected by electro-reflectance measurements for the sam- 
ples with nitrogen concentrations of >0.8% 0,0 ■ This 
peak is located at about 0.4-0.8 eV above the conduction 
band (CB) minimum (denoted as EL) and shifts upward 
in energy with increase of nitrogen concentration [3j,|4j,|5|. 

Several physical mechanisms have been proposed to 
describe these nitrogen-induced changes : two-level band 
anticrossing (BAC) model |6j,l7J, disorder-allowed T-X- 
L cou plin g; [g, 0, and the formation of the impurity 
band 10] . In the frames of the BAC model the E+ 
state naturally appears as the high-energy solution of 
a two-state Hamiltonian, describing the interaction of 
the localized nitrogen state with the extended CB states 
of the host matrix compound [Til ]. However, the first- 
principles EI E3 as well as empirical pseudopo- 
tential [illlialla] calculations failed to validate the BAC 
model. On the other hand, based on these calculations, 
the intraband coupling model has been proposed. In this 
model the E + peak is caused by a disorder-allowed tran- 
sition from the valence band (VB) maximum to the L- 
point of the conduction band. In particular, it was pro- 
posed that the energy position of the E + transition is 
a configuration-weighted average of a nitrogen impurity 
state and CB L-state, denoted as L\ c 0. However, to our 
knowledge, there are no studies of the GaAsN system, 
which identify and clarify the origin of the E + transition 
based on direct first-principles calculations of the optical 
properties of this system. 



In this Letter we report the results of theoretical study 
of the optical properties of GaAsN system using state- 
of-the-art ab-initio method based on density functional 
theory (DFT). The E + transition is identified by direct 
calculation of the dielectric function of the system for 
nitrogen concentrations of 0.9-3.7%. Its origin is further 
clarified by studying band-to-band contributions to the 
dielectric function. To understand the effects of atomic 
ordering, the supercells of different symmetry are con- 
structed and studied separately. Our results show that 
the E + transition is formed by the interplay of two contri- 
butions : the T-point transitions to the L lc -derived states 
of the GaAs matrix, and the L-point transitions to the 
impurity-modified first conduction band of the system. 
Further, we study the changes in these contributions, 
caused by different impurity distribution, and show the 
importance of nitrogen potential in the formation of the 
E + transition. 

Our calculations were performed within the local den- 
sity approximation (LDA) and generalized gradient ap- 
proximation (GGA) to the density functional theory |17| . 
Four supercells were constructed to describe GaAsi-zNr 
system with nitrogen concentrations in the range of 
x =0.0093-0.037. In 64-atom and 216-atom supercells ni- 
trogen atoms form a simple cubic (SC) lattice, while in 
54-atom and 128-atom supercells a face-centered cubic 
(FCC) impurity sublattice is formed. A pseudopotential 
planewave method as implemented in ABINIT |l8| code 
was employed for full geometry optimization of the su- 
percells. This included structural relaxation both with 
respect to the volume of the unit cell and to the in- 
ternal coordinates of the atoms. A 16 Ry basis cutoff 
and 2x2x2 Monkhorst and Pack ^1 sampling of the 
Brillouin zone (BZ) showed good results for structural 
relaxation. The LDA exchange-correlation functional of 
Perdew and Zunger |20j has been used for calculations. 

Having obtained geometrical characteristics of the su- 
percells, we proceeded further with a highly-accurate lin- 
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Fig. 1 - Calculated imaginary part of the dielectric 
function for GaAsi^xN^ systems, normalized per unit 
volume. A progressive Lorentzian smearing with full 
width at half maximum (FWHM) of 0.01(?iw) 2 has been 
applied. 

earized augmented plane wave method (LAPW) as imple- 
mented in WIEN2K program package |2j|. This method 
belongs to the family of all-electron DFT-type methods 
and does not use any pseudopotentials. Band structures 
and optical properties of the model systems have been 
calculated using the GGA exchange-correlation func- 
tional of Engel and Vosko |2^. This functional signifi- 
cantly improves the values of the band gap and effective 
electron and hole masses for GaAs |23j as well as for 
GaAsN [24| systems. For calculation of optical proper- 
ties a dense 20 x 20 x 20 k-grid has been used [25( . Having 
obtained selfconsistently the crystal wavefunctions |kn) 
and band energies E^ n for n bands of the system, the 
imaginary part of the dielectric function £2(w) can be 
calculated by summing direct transitions over the BZ, 
weighted with corresponding matrix element of the tran- 
sition probability : 

nn' 

X/kn(l - fan>)5(Ek n > - E kn - Hoj). 

Here p=—ihV is the momentum operator, and /kn is 
the Fermi distribution function, which ensures that only 
transitions from occupied to unoccupied states are in- 
cluded in summation. 

We present in FigQ calculated imaginary part of the 
dielectric function (£2) for GaAsN systems with nitro- 
gen concentrations of 0.93, 1.56, 3.13, and 3.70%. A cal- 
culated £2 for pure GaAs is also presented for sake of 
comparison. All main features, observed experimentally 
for GaAs and GaAsN systems, are well-reproduced by 
theoretical calculations. The common features of the di- 
electric function for the pure and doped systems include 
the presence of the fundamental absorption edge Eq, as 
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Fig. 2 - Optical transition energies vs band gap 
reduction : the Eq, E + , and E\ transitions are 
presented. 

well as two main maxima, namely peaks E\ and E-i- The 
origin of these features is well-understood and described 
in the literature (2|| . The fundamental absorption edge is 
due to transitions between valence band maximum and 
conduction band minimum, which occur in the center 
of the Brillouin zone. The E\ peak is related to direct 
transitions along the eight equivalent T — L directions of 
the BZ. And, finally, the absolute maximum E2 contains 
contributions from the direct transitions occurring near 
A-point of the BZ. 

It is clearly seen that with increase of nitrogen dop- 
ing the fundamental edge moves down to lower energies. 
This reflects the effect of the nitrogen-caused band gap 
reduction, and agrees well with experimental observa- 
tions 0, I3| • We also notice that the E 2 peak hardly 
moves with increase of the nitrogen contents, while the 
Ei transition moves slightly to higher energies. The main 
feature, which makes the spectra of GaAs and GaAsN 
systems substantially different, is the appearance of a 
new nitrogen-induced peak, which we attribute to the 
E + transition, observed experimentally With increase of 
nitrogen content the E+ peak shifts to higher energies 
faster than the E\ transition. 

To represent our results quantitatively, we plot in FigEl 
the optical transition energies as a function of band gap 
reduction i?o(GaAs)-i?o (GaAsN) with the energies mea- 
sured relative to the VB maximum. In this representation 
the Eq transition is a straight line with a slope of -1.0. 
The energy evolution of the E\ peak can also be well- 
approximated by a straight line with a slope of +0.25. 
This slope is close to +0.22, experimentally obtained by 
Shan etal [ilj for four different nitrogen concentrations 
with further parametrization of the results within the 
BAC model. 

While for Eq, Ei and E2 the exact transition energies 
can be easily determined directly from the dielectric func- 
tion, this task becomes much more difficult for the E + 
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Fig. 3 - (a) and (c) : Contributions to the imaginary 
part of the dielectric function from different band 
combinations, calculated for 128- and 64- atom 
supercells. (b) and (d) : Calculated electronic 
bandstructure of these supercells showing the various 
direct transitions, responsible for the features in £2. 



peak. This transition is clearly visible only for 64- and 
128- atom supercells. For the system with the highest ni- 
trogen concentration of 3.7% this transition is smeared 
to large extend probably due to strong N-N interactions, 
while for the smallest concentration of 0.93% it is hardly 
visible due to a reduced oscillator strength. By comparing 
the dielectric function for 128- and 64- atom supercells 
with corresponding energy spectra at the T-point, we find 
that the E+ transition energy for these structures exactly 
corresponds to the position of the a\(Li c ) singlet. At this 
point we are making an important assumption, which will 
be later verified by a detailed analysis. We assume that 
for the other two structures the E\ peak also originates 
from a T transition to the a\{Li c ) singlet. This assump- 
tion is also supported by the fact that for our model 
system, containing 0.93% of nitrogen, this singlet is lo- 
cated at 0.53 eV above the conduction band minimum, 
which is in good agreement with experimental value of 
0.57 eV, obtained for the E+ peak for the samples with 
1% of nitrogen 0- 

Now we plot in Fig|21the energies of the E+ transitions 
within the same representation. Linear approximations, 
performed separately for SC- and FCC- types of super- 
cells, showed practically the same slope of +0.40, which 
is smaller than experimental value of +0.67, obtained by 
Perkins et al [j| for the samples with four different ni- 
trogen concentrations. Extrapolation to the dilute limit 
gives the average value of E + — Eq(x = 0) = 145 meV, 
which is close to the experimental interval of 150-180 
meV, obtained for the energy position of nitrogen level in 
GaAs 0| . However, such a detailed comparison should 
be done with a caution due to the fact that only two 
points are available for each type of structures. 



We now try to understand the origin of the E + tran- 
sition in the GaAsi_ x Na; system. The expression, used 
for the calculation of £2(0;), gives us a natural possibility 
to calculate the contributions to the dielectric function 
from diffrent band combinations. FigureEJa,c) represents 
this decomposition of £2, performed for the systems with 
1.56% and 3.13% of nitrogen, taking into consideration 
only three top valence and five lowest conduction bands. 
Only these combinations were found to be responsible 
for the formation of the dielectric function in the energy 
interval 0.0-2.5 eV, which includes the E+ transition. 

Let us first concentrate on the system with 1.56% of 
nitrogen (Fig. Et). We consider three contributions to 
the dielectric function : transitions to the first conduc- 
tion band (Ci), to the second conduction band (C2), and 
cumulative transitions to the conduction bands 3, 4 and 
5 (C345). We observe that the transitions to C345 do not 
contribute to the near-edge features of £2. The results 
show that the E + feature is formed by contributions of 
two transitions, marked A and B, which both are energet- 
ically located near 1.65 eV, but have completely different 
origin. Analyzing the energy positions of peaks A and B, 
and keeping in mind that they come from the transitions 
to different conduction bands Ci and C2, we are able 
to uniquely connect these transitions with corresponding 
saddle points in the bandstructure of the system. Figure 
Et> represents the bandstructure of the 128-atom super- 
cell, showing three top valence and five lowest conduction 
bands. The results show that the peak A is formed by the 
transitions from the the valence bands to the first con- 
duction band near the L-point of the BZ, while the peak 
B is due to the transitions to the second conduction band 
near the T-point. The fact that these two transitions are 
very close in energy and have significant intensities, leads 
to the formation of a well-defined feature in the spec- 
trum, which we attribute to be the E + transition. The 
transitions near the A"-point of the BZ contribute to the 
spectrum in the form of the peak C. It is formed by the 
transitions to the second conduction band, but it is twice 
smaller than peak B in intensity, and is located at higher 
energies. 

It was found that large configuration-induced bandgap 
fluctuations are present in the GaAsN systems both in 
cases of one-dimensional nitrogen chains 1281. as well as 
at relatively low nitrogen concentrations 24]. It seems 
interesting to study the influence of impurity ordering 
on the formation of the E + transition. So, we proceed 
further with the same decomposition analysis of £2 for 
the 64-atom supercell, where the nitrogen atoms form 
a simple cubic sublattice. The contributions of the Ci, 
C2 and C345 bands are represented in Fig. [3];, and the 
band structure with corresponding optical transitions is 
shown in Fig. EH- The results show than the E+ transition 
is also formed by a combination of peaks A and B, but 
the relative intensities of these peaks change. Comparing 
Fig. Et an d Fig- Et we observe that the contribution of 
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peak B decreases, while the one of peak A increases. The 
increase of the contribution of the band Ci for the SC- 
type supercell can be explained by the appearance of the 
additional saddle point in the T — R valley. The saddle 
point of the same type appeares also in the C2 band, 
which leads to the formation of additional peak D, which 
is not observed in the FCC-type structure. 

We now conclude our study by exploring the role of im- 
purity potential in the process of the formation of the E+ 
transition. Besides splitting the instates into a singlet 
and a triplet, the introduction of nitrogen also leads to 
a strong perturbation of the first conduction band of the 
host compound 0|. These two effects could be caused 
by both the geometrical deformations near the impurity 
atom, and by the nitrogen potential itself. In this last 
part of our study we are making an attempt to verify the 
influence of these two effects on the formation of the E + 
transition. 

To gain further insight we performed calculations of 
the electronic structure for two model systems with sub- 
sequent analysis of different contributions to the dielec- 
tric function. We have first taken the 128-atom supercell 
with the fully relaxed geometry and replaced nitrogen 
atom back to arsenic. Thus, in this model we have re- 
moved the effects of the impurity potential, staying only 
with the effects of structural relaxation. The results show 
two dramatical changes in the electronic structure and 
optical properties of the system : the band gap of this 
model structure opens to the value of 1.29 eV, which is 
close to 1.26 eV, obtained for the GaAs structure, and 
the E + transition is no longer present in the spectrum. 
The contribution analysis demonstrated the reduction of 
the amplitudes of peaks A and B to practically negligible 
values. Moreover, in spite of local geometrical deforma- 
tions, the splitting of L\ c states is now only 0.06 eV, 
which is almost twice lower than 0.11 eV, observed for 
the N-doped supercell. This demonstrates that local me- 
chanical deformations introduce only small perturbations 
to the GaAs system, which are not responsible neither for 
dramatic band gap reduction, nor for the appearance of 
the E+ transition. 

These conclusions are supported by our results, ob- 
tained for the second model structure. In this case we 
have taken the N-doped 128-atom supercell, but stayed 
within the geometry of the GaAs structure. This model 
allows us to observe only the effects of the nitrogen po- 
tential without any influence of impurity-caused defor- 
mations. The decomposition analysis showed that both 
peaks A and B, which are responsible for the formation 
of the E + transition, are present in the spectrum. The 
differences from the relaxed structure include the larger 
splitting of the L\ c states (0.17 eV) and shifts in the en- 
ergy positions of peaks A and B. We also observe open- 
ing of the band gap up to 1.13 eV. However, this value 
is smaller than the one, obtained for the previous model 
system, where only the absense of the nitrogen potential 



caused a dramatic opening of the band gap to the one of 
the GaAs structure. 

In conclusion, we have presented theoretical study of 
the optical properties of the GaAsN system. The origin 
of the experimentally observed E+ transition is identified 
by the analysis of the dielectric function, calculated ab- 
initio for four different supercells. Our results show that 
the main contribution to the E+ peak is given by T-point 
transitions from the VB maximum to the ai(Li c )-singlet, 
which originates from a nitrogen-induced splitting of the 
CB L-point. At the same time, nitrogen- induced modi- 
fications of the first CB are found to give rise to addi- 
tional L-point direct transitions. Being close in energy 
to T-point transitions, they are shown to contribute sig- 
nificantly to the formation of the E + peak. Finaly, local 
geometrical deformations are found to be only a second 
order effect in the process of the formation of the E + 
transition, and the influence of the impurity potential it- 
self is shown to be responsible for the E + transition and 
dramatic band gap reduction in GaAsN. 
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